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We construct, by numerical means, static solutions of the spheri- 
cally symmetric Einstein- Vlasov-Maxwell system and investigate vari- 
ous features of the solutions. This extends a previous investigation [5] 
of the chargeless case. We study the possible shapes of the energy den- 
sity profile as a function of the area radius when the electric charge of 
an individual particle is varied as a parameter. We find profiles which 
are multi-peaked, where the peaks are separated either by vacuum or 
a thin atmosphere, and we find that for a sufficiently large charge pa- 
rameter there are no physically meaningful solutions. Furthermore, we 
investigate if the inequality 



derived in [5], is sharp within the class of solutions to the Einstein- 
Vlasov-Maxwell system. Here M is the ADM mass, Q the charge, 
and R the area radius of the boundary of the static object. We find 
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two classes of solutions with this property, while there is only one 
in the chargeless case. In particular we find numerical evidence for 
the existence of arbitrarily thin shell solutions to the Einstein- Vlasov- 
Maxwell system. Finally, we consider one parameter families of steady 
states, and we find spirals in the mass-radius diagram for all examples 
of the microscopic equation of state which we consider. 

1 Introduction 

In this work matter is described as a large ensemble of charged particles 
which interact via the gravitational and electromagnetic fields created by 
the particles themselves. All the particles have the same rest mass, normal- 
ized to 1, and the same charge qo > 0. The distribution of the particles on 
phase space is given by a density function /. The particle ensemble is as- 
sumed to be collisionless which implies that / satisfies the Vlasov equation. 
Macroscopic quantities such as mass-energy density, pressure, and charge 
current, which act as source terms in the field equations, are obtained by 
integrating / with respect to specific weight functions. The resulting system 
is called the Einstein- Vlasov-Maxwell system and is stated in Section [21 In 
the present work we construct, by numerical means, static solutions of the 
asymptotically flat, spherically symmetric Einstein- Vlasov-Maxwell system, 
and we investigate three different features of the solutions. In the chargeless 
case, i.e., for the spherically symmetric Einstein- Vlasov system, a similar 
study has been carried out in [5]. There is an essential difference between 
these two systems concerning known mathematical results of existence of 
static solutions and their properties. In the chargeless case it is known that 
a wide variety of static solutions with finite extent and finite ADM mass 
exist, cf. [U [H], and the references therein, whereas the problem of exis- 
tence of static solutions in the charged case has not yet been studied. The 
mathematical construction of steady states in the chargeless case is based 
on a certain ansatz for the density function /. Here this ansatz is modified 
to handle the charged situation, cf. Section O We believe that this con- 
stitutes a natural starting point for showing existence of static solutions of 
the Einstein- Vlasov-Maxwell system, but we do not include such an analy- 
sis here since the purpose of the present paper is to investigate numerically 
three features of static solutions which we now describe in some detail. 

In Section 4 an analysis of the behavior of the energy density as a function 
of area radius is carried out for different values of the charge parameter qq. 
Qualitatively we find a similar structure as in the chargeless case [5J, e.g. 
there are solutions with an arbitrary number of peaks, and these peaks are 
separated either by vacuum or by a thin atmosphere. The choice of charge 
parameter affects in some cases the number of peaks. If the charge parameter 
reaches a certain critical value the solutions break down before the energy 
density vanishes. It is natural to compare this with the Newtonian Vlasov- 
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Poisson system for which there is no difference in the form of the equations 
whether one models a mono-charged plasma or a gravitating system, except 
for the sign in front of the force field. If one studies a charged gravitating 
system this sign is positive or negative depending on whether go < 1 or 
go > 1- Steady states with finite extent only exists in the former case where 
the effective force field is attractive. This is in accordance with what we find 
in the relativistic situation. 

In Section 5 we investigate if there are static solutions of the Einstein- 
Vlasov-Maxwell system such that the inequality 



which was derived in [2J, can be saturated in the sense that the quotient of 
the left and right hand side in (jl.ip is arbitrarily close to one. If there are 
such solutions we say that the inequality (jl.ip is sharp within this class of 
solutions. Here M is the ADM mass, Q the total charge, and R the area 
radius of the boundary of the static object. It was shown in [2j that (jl.ip 
holds for any static solution of the Einstein-Maxwell-matter system which 
satisfies the energy condition 



where p > and ^ are the radial and tangential pressures respectively 
and yO > is the energy density. This condition is satisfied by Vlasov matter. 
A suitable extension of the inequality (jl.ip also holds inside the static object, 
cf. ^ and Section 5. Moreover, it was shown in [2] that the inequality is 
sharp, and in particular that equality is attained by infinitely thin shell 
solutions. The method of proof in [2j is quite general and applies to any 
matter model for which (jl.2p holds. However, the solution constructed in 
the analysis leading to sharpness has features which solutions of the Einstein- 
Vlasov-Maxwell system do not have. Hence, one motivation for the present 
study is to investigate if sharpness of the inequality can be attained by 
solutions when a real matter model is chosen so that the system of equations 
includes a matter field equation, in the case at hand the Vlasov equation. 
In contrast to the chargeless case it is not known if there are solutions 
other than infinitely thin shells which saturate (jl.ip . and it is not known 
if arbitrarily thin shell solutions do exist for the Einstein- Vlasov-Maxwell 
system. 

In the uncharged case more is known. It follows from [3] that infinitely 
thin shell solutions are unique in saturating the inequality (jl.ip for (5 = 0, 
i.e., the inequality 

4ii 

M < — . 
- 9 

Since an infinitely thin shell solution is not a regular solution of the Einstein- 
matter system this statement should be interpreted in the sense that a se- 




(1.1) 



P + PT < P, 



(1.2) 
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quence of regular solutions tending to an infinitely thin shell will in the 
limit give equality. Moreover, it is known [3] that regular arbitrarily thin 
shell solutions of the Einstein- Vlasov system do exist, which then in partic- 
ular implies that there are steady states to this system such that M/R is 
arbitrarily close to 4/9. 

In the present work we find numerical evidence for answering the issues 
raised above. Indeed, we construct arbitrarily thin shell solutions to the 
Einstein- Vlasov-Maxwell system, which saturate the inequality (jl.ip in the 
limit. Moreover, in contrast to the uncharged case we also find another 
type of solutions which saturate the inequality. These solutions have the 
feature that M, Q, and R are all equal; they represent an extremal object. 
The latter property may be of interest in fundamental black hole physics, 

cf. [giinii. 

In Section 6 the third and final property is investigated, namely the 
relation between the ADM mass and the outer area radius of a one parameter 
family of steady states to the Einstein- Vlasov-Maxwell system. The one 
parameter family is obtained by prescribing the way in which / depends on 
the local energy and the angular momentum, which we call the microscopic 
equation of state. We find numerical support for mass-radius spirals for all 
examples of the microscopic equation of state which we investigate. This 
agrees with the result in [5] in the chargeless case, but on the other hand 
it differs from the Newtonian situation where the presence of such spirals 
heavily depends on the microscopic equation of state. In [HI [13] the question 
of which equations of state in the fiuid case give rise to spirals is investigated. 

To conclude this introduction we mention [6l [HI [10] where numerical 
investigations of solutions to the Einstein-Maxwell system are carried out 
using other conditions on the matter model. These studies focus on the issue 
of bounding M in terms of Q and R. 

2 The spherically symmetric Einstein-Maxwell- 
Vlasov system 

We choose general local coordinates x° on the spacetime manifold and we 
denote by the corresponding canonical momenta; Greek indices always 
run from to 3 and Latin ones from 1 to 3. We assume that = t \s a. 
timelike coordinate and that can be expressed by through the condition 
that all the particles have rest mass normalized to 1: gai3P°'P^ = — 1- Then 
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the Einstein- Vlasov-Maxwell system takes the following form: 

F°^.„=47rJ^, (2.1) 

^/37;a + -^7q;/3 + -^q:/3;7 = 0) (2-2) 
1 

Ra(3 - -j^aapR = 8TT{Ta/3 + Tq/j), (2.3) 

dtf + ^d,.f - ^(ri^pV + qop''Fj)d^.f = 0, (2.4) 



where 



f — d?p 

Ta/3 = - PaPfsVlal , (2.5) 

rap = ^ {F^'Fp, - ^F,,F^^) , (2.6) 

j' = ciJ //v¥^. (2.7) 



Here ()2.ip and ()2.2p are the Maxwell equations, ()2.3p are the Einstein equa- 
tions, (|2.4p is the Vlasov equation, and F-a denotes the covariant derivative. 

We consider this system under the assumption of spherical symmetry. 
Hence the metric, expressed in Schwarzschild coordinates, takes the form 

where 

t > 0, r > 0, 6* G [0, vr], G [0, 27r]. (2.9) 

For the metric to approach that of Minkowski space as r goes to infinity, 
the boundary conditions 

lim \{t,r) = lim /x(t,r) = (2.10) 
r— >oo r— >oo 

are imposed. Furthermore, the condition 

A(t,0) = (2.11) 

ensures a regular center. We introduce the corresponding Cartesian coordi- 
nates X = {x^,x^,x'^) = r(sin 0COS (/?, sin6' sin(^, cos 6*) and find that 



Po = -e^A/l + |p|2 + (e2A-l) ^ 



.2A .m'^-P 



Here x ■ p denotes the Euclidean scalar product of the vectors x = 
(x^, x^, rr'^), p = {p^,p'^,p^), and | • | denotes the Euclidean norm on M^. 
For a spherically symmetric electric field of the form 



E' = e- 
r 
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the non-zero components of the electromagnetic field-strength tensor are 



?0i 



X 



r 



Let 



where 



X ' V . .0 2/1 i2 1\ 

w = , L = \x X v\ =r{\v\ —w), 

= + {e^ - ly- 



r r 

the variables w and L can be viewed as the momentum in the radial direc- 
tion and the square of the angular momentum, respectively, expressed in a 
suitable frame. The system now reads 

q' = Anr^pg, (2.12) 

e-2^(2rA' -!) + ! = Swr^p + \, (2.13) 

e~^^(2r//' + 1) - 1 = Sirr'^p - ^, (2.14) 



dtf + e^-^ ^ — drf - f e^^- Wl + + L/r^ + Xw 



_eM^ - e'^-^ , ] d^f = 0, (2.15) 



where 



roo poo 

p=-/ / ^l + w-^ + iy/dLdw, (2.16) 

1^ J-oo Jo 



-oo JO 



P=^/ / -j====fdLdw, (2.17) 



oo POO 



Pq = qoe^^ I I fdLdw. (2.18) 

^ J-oo Jo 

Here a prime or dot denotes the derivative with respect to r or t, respectively, 
Pg = Pq{t-,^) is the charge density, q = q{t,r) is the charge contained in the 
ball with area radius r about the origin, p = p{t, r) is the energy density as 
defined when no charge is present, p = p(t, r) is the radial pressure, and the 
modulus of the electric field is given by 
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3 Constructing static solutions to the spherically 
symmetric Einstein-Vlasov-Maxwell system 

In this paper we are interested in static solutions, so ()2.15p reduces to 

'^^rf+ (^e^^^/l + w^+L/r^ + ^ - ^'(1 + + L/r^)^ d^f = 0, (3.1) 

where / = f{r,w,L), fi = fj,{r) and A = A(r). Due to spherical symmetry 
the quantity L is conserved along characteristics of the Vlasov equation, and 
so is the particle energy E defined as 

E = e'^Vl + u'2 + L/r2- go F e2^(^)+''(^)e(7?)(i7/. (3.2) 

^0 

Hence any density function of the form 

f{r,w,L) = '^{E,L) (3.3) 



satisfies the static Vlasov equation ()3.ip . In order to motivate (13. 2p we 
combine the electric potential (pE and the magnetic vector potential A into 
a four-vector 

K° = = A\ 

The electromagnetic field-strength tensor can be derived as 

In particular, 

For = -e^^+^e = -drKo, 
i.e., with the electric potential taken to be zero at r = 0, 

KQ = re2^(^)+'^(^)e(7?)dr/. 

^0 

From this we get the particle energy E = — (po + 90/^0) as defined in ()3.2p . 

The ansatz (j3.3p is a generalization to the charged case of the standard 
ansatz for the Einstein- Vlasov system which is obtained for go = 0, and to 
the best of our knowledge it has not appeared in the literature. Equations 
()2.12p -- (l2.14ll can be rewritten as the system of ODE's 

'^-q = Anr'^pq (3.4) 



dr 

re"2^) = 1 - SvrrV - ^ (3.5) 

(3.6) 



dr \ / 
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where the quantities Pq,p^ and p are now functionals of g,A, and p. In 
order to obtain a steady state with finite ADM mass and finite extension 
we prescribe some cut-off energy Eq > {) and assume that ^{E,L) = for 
E > Eq. Taking this into account, 



2yj- /■'Umax rLmn^ 

p=—j / v^l + ^2 + L/r2$ dL dw, 



/o Jo 

p=^ / ^ dL dw, 

A 



Pq = qoe — / ^dLdw, 



JO 



1 

2 \ 2 



where the upper hmits 

w^max = (^6-2/^ (^i?o + 'Zo e^^^^^+^^^h{r^)dri^ - 1 j 
and 

follow from the condition E < Eq. Since e(0) = A(0) = 0, (l33I)-(lM|) can be 
solved if p{0), qo, and Eq are specified. In order to continue we introduce 



2/1 = e^e = — , y2 = e^^, ys = (3.7) 
r tjQ 



By making the ansatz ^{E,L) = (f){E/Eo,L), ([Till - lfT^ turn into 



,r2 

d f r 



r~yi) = 4nr^Pg, (3.8) 



) = 1 - (8vr/5 + yj) , 



, , , ^ . v-r . , (3.9) 
^ {rys) = 2/2^3 (1 + (8vrp - y?)) , (3.10) 



where 



2tt 



r^max /'-^'max 

/ / yrTw2TX7^</'dLdu;, (3.11) 
Jo ^0 

p=^ / (pdLdw, (3.12) 

^'70 Jo + + 

Pq = qoVy^— / (l)dLdw, (3.13) 
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with upper limits for the integrals given by 



Thus, since 



En 



V r~ jQ 



i + w"^ + ^ - qo / Vy2iv)y3iv)yii'n)dr] 



the dependence on £"0 is eliminated from the system of ODE's (|3.8|) - (|3.1U|) . 
and the equations can be solved after specifying merely 2/3(0) and go- Notice 
that from the boundary conditions (|2.1U|) we require that e2A'('^) ^ 1 as 
r ^ 00 so that the limit of 2/3 at infinity and (13. 7p determine Eq. 

For almost all numerical solutions presented below an ansatz of the form 

0(ry,L) = (l-r?)^(L-Lo)^, (3.14) 

where k > 0, I > —1/2, k < 31 + 7/2, Lq > 0, has been used. Here 
X-i- = maxja;, 0}. Our code does allow for a wide variety of ansatz functions, 
but the qualitative behavior is similar in all cases we have tried, cf. the end 
of Section 6. When using an ansatz with a cut-off Lq > for the square of 
the angular momentum as in (j3.14p . then at any point with 



r < 



\ 



^ (i + qo fo y2{ri)yz{n)yi{n)dn 



\ 



the matter quantities ()3.1ip ~ ()3.13p are zero. In particular the matter quan- 
tities will be zero and the functions yi,y2iy3 will be constant for r < Rq 
with 



Ro 



Ln 



y3(o)-i-i" 



The integrals in the matter quantities are calculated using the piecewise 
Simpson's rule, and the differential equations are solved, from r = Rq radi- 
ally outwards, using the Euler method with a variable step length /i„, given 
at a point r„ as 

hr. 



ln(e + |p^|)' 

where p'^ denotes the second derivative with respect to r at r = r„, /imax 
is an appropriately chosen maximum step length, and e is Euler's number. 



9 



This implies that the code resolves those regions more finely where p varies 
rapidly. In the chargeless case the integrals in (j3.1ip and (j3.12p can be 
carried out explicitly for suitable ansatz functions cf. [5, and we used 
the corresponding code to test the present one where these integrals are 
evaluated numerically. 

4 Characterization of steady states 

In this section we study the behavior of the uncharged energy density p as a 
function of the area radius when the parameter for electric charge is varied. 
We find profiles which are multi-peaked, where the peaks are separated either 
by vacuum or a thin atmosphere. In this respect our results are qualitatively 
similar to the chargeless case, but the value of the charge parameter does 
affect the number of peaks in some cases. On the other hand, we find that 
for a sufficiently large value of the charge parameter there are no physically 
meaningful solutions. This is intuitively clear in view of the Newtonian 
case where no steady states of finite mass exist when the effective force is 
repulsive. 

The steady states for charged matter largely follow the same basic struc- 
tures as noted in [5] for the non-charged case, i.e., we get steady states with 
support for p in [Rq, i?], i?o > 0, called shell configurations, and states where 
Rq = 0, called ball configurations. 

To visualize the steady states, for each quadruple {k, Lqi 2/3(0)) a figure 
containing two subfigures is presented. Subfigure (a) shows starting points, 
stopping points, i.e., the inner and outer boundaries of the matter shells, as 
well as maxima of p for e [0, (?o,c[) while subfigure (b) shows individual 
solutions for three values of q^, namely go = 0,0.5, 1.0. The starting points 
in subfigure (a) are shown as dashed lines, the stopping points as solid lines 
and peaks are represented by dotted lines. The plots complement each other, 
the plot for the starting and stopping points contains little information on 
the shape of the steady states, while vacuum regions can be difficult to 
notice in the plots for individual solutions. To remedy the fact that peaks 
for p closer to r = are in general much larger in magnitude than peaks 
farther out, ln(47rr^/? + 1) rather than p has been plotted against r for the 
individual solutions. It should be borne in mind that this has the effect that 
the positions of (and in one case also number of) maxima in subfigures (a) 
and (b) differ. 

The most noticeable effect of changing the particle charge q^ is that as 
go approaches a critical value go,c, the outer radius R of the support for p 
increases dramatically, as can be seen in Figures [T] and Intuitively, this 
is expected, since for some value of qo, the repulsive forces between individual 
particles from electric charge balances the attractive forces of gravity. For 
values of go larger than go,c) the numerical solution breaks down since at 
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q^=1, k=0, 1=1.5, L^=0.2, y3(0)=0.05 

1.4 r 



0.2 - 

- 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 



Figure 2: In (47rr^p + 1) and In (47rr^pq + 1) for a shell configuration 

some point y2{'f) approaches infinity. Thus we cannot have arbitrarily large 
particle charge and still obtain a solution. In Figure [2] the behaviour of 
In {Anr'^p + 1) and In (47rr^pg + 1) is displayed. Recall here the definition of 
Pq in p.lSp . As can be seen the profiles of these quantities are similar. This 
observation applies to other cases in this paper as well which is the reason 
why we only focus on In (Airr'^p + 1). 

The value of go,c for fixed values of k,l and Lq varies with 1/3(0) as can 
be seen in Figure El We see that qo^ciusiO)) has an undulating quality with 
increasing amplitude and decreasing frequency of oscillation for larger 1/3(0). 
As 1/3(0) 1, we do always have that go,c — > 1- This is easily understood, 
since in the Newtonian case the critical particle charge is exactly 1, and 
in the limit 1/3(0) — > 1 the solutions become essentially Newtonian; notice 
that 1/1/3(0) — 1 is the central redshift factor [SJ Eqn. (2.20)] which is a 
measure for how relativistic the system is. At any point where p equals 
zero, a Reissner-Nordstrom solution with appropriate initial values can be 
joined to form a steady state with finite support. The above limit for qq 
thus only applies to situations where the distribution function is given by a 
single expression on the entire domain interval r G [0, 00 [. 

Figure [T] shows a double-peaked, single-shelled shell configuration (k = 
0,1 = 1.5, Lq = 0.2,1/3(0) = 0.05). The positions and magnitudes of peaks 
are virtually unaffected by the value of qo and the only effect that can be 
seen is that the tail grows in length and magnitude. The outer radius R 
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Figure 4: Shell configurations, decreasing R 
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y3(0) 



Figure 6: go,c as a function of ^3(0) 

increases strictly monotonically. In Figure [3] a double-peaked ball configu- 
ration with parameters (1/3(0) = 0.01, /c = 1,/ = 12, Lq = 0) is displayed. 
Here we see that as qo increases, an extra maximum of p appears. In this 
case R does not increase strictly monotonically, as a decrease can be seen 
before the outer radius finally blows up as qo approaches go,c- Solutions in 
which a new shell appears and with strictly increasing R can however be 
constructed. Although barely noticable in Figure [3l the radial position of 
the new maximum that appears for larger qQ will increase at first and then 
decrease. This effect is more pronounced in Figure HI however. The afore- 
mensioned case is displaying the same behavior as in Figure [3l this time for 
a single-peaked shell configuration (1/3(0) = 0.16, A; = 0,1 = 1.5, Lq = 0.2). 
For higher values of Lq multi-shelled configurations (i.e., configurations with 
multiple peaks, separated by vacuum regions), can be obtained. The effect 
on these is that as qo increases, one or more additional shells appear as can 
be seen in Figure E] (2/3(0) = 0.0025, /c = 1,/ = 3,Lo = 100). These newly 
appearing shells mimic the behavior of the newly appearing peak in Fig- 
ure HI In Figure [5] it can be seen clearly that all peaks, except the innermost 
one, that are present at go = are showing an inclination to move towards 
r = 0. This behavior is also present in all cases with more than one peak, 
although not noticable in the plots. The innermost peak, on the other hand, 
has a contrary tendency to move outwards. 
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5 Sharpness issues of the main inequaUty 



The purpose of this section is to investigate aspects concerning sharpness of 
the inequahty 



where rUg is the total gravitational mass given by 

mg{r)=47rj rf p{ri)drj + {j -^-^drj + j =: rui + niq. (5.2) 

This inequality was derived in [2j and it was shown to hold for any static 
solution of the Einstein-Maxwell- matter system which satisfies p.2p . In 
addition it was assumed that the solutions satisfy 



q < rUg, nig + \ rn^, — <r. (5.3) 



The latter conditions are imposed to ensure that the solutions are physically 
meaningful, cf. [10] . To better understand the motivation for our study we 
recall the results in the uncharged case. 

If Q = the inequality ()5.ip reduces to the Buchdahl inequality 

^ S \, (5.4) 

which was first proved in |7J under the Buchdahl assumptions that the pres- 
sure is isotropic and the energy density is non-increasing outwards. The 
inequality was then shown to hold independently of the Buchdahl assump- 
tions [3] for solutions which satisfy the energy condition p + 2pT < p. A 
different proof was later given in [12]. The advantage of the latter method 
is that the proof is shorter and more flexible since it allows for other energy 
conditions than (II. 2p . The disadvantage lies in the issues of sharpness and 
construction of the saturating solution. Firstly, the method does not imply 
that the class of saturating solutions is unique, and secondly, it is not clear 
that a solution to a coupled Einstein-matter system can have the proper- 
ties of the saturating solution constructed in [12]. In particular solutions 
to the Einstein- Vlasov system are ruled out. These issues have however 
affirmative answers. Uniqueness is obtained in [3] where it is proved that 
a saturating solution must be an infinitely thin shell solution. In [4J it is 
shown that regular, arbitrarily thin shell solutions of the Einstein- Vlasov 
system exist, which implies that there are steady states to this system with 
M/R arbitrarily close to 4/9. 

Let us now return to the charged case and discuss the latter two issues. 
The proof of (|5.ip in [2] is based on the method in [T2] and the method 
in [3] does not apply. A proof based on the latter method would imply 
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uniqueness of the saturating solution, and it is thus natural to ask if there 
is another class than infinitely thin shell solutions with this property. Also, 
although the result in [2] shows that infinitely thin shell solutions saturate 
the inequality there is no analogous result to [3] in the charged case, i.e., 
the question whether or not the Einstein- Vlasov-Maxwell system admits 
arbitrarily thin shell solutions is open. Below we will present numerical 
evidence that the system does admit arbitrary thin shell solutions and in 
addition that another class of saturating solutions does exist. 
We introduce the quantity 

y/mg{r) 
r := sup — ===, 

r>o y/r r q [r) 

which in terms of the functions introduced in ()3.7p reads 

3^2 \/ 1 - y2Hr) + 47rr2y2(r) 

r = sup , — . 

r>o 2 1 + y/l + 127rr2y2(r) 

By (jS.ip . r is subject to the inequality 

r < 1. (5.5) 

In Figure [7] the numerical results for a shell configuration with (k = 1,1 = 
1, Lq = 10) are displayed as follows. F as a function of 2/3(0) is displayed in 
ascending order for 

go G {0, 0.25, 0.5, 0.75, 0.999, 1.001, 1.25, 1.5}, 

the bottom-most curve being that for qq = 0, the next to bottom-most being 
that for Qo = 0.25 and so on. In all cases tested for other parameter values, 
r (7/3(0)) falls within the first shell of the solution. 

We see that as 2/3(0) approaches zero (15. 5|) approaches equality for all val- 
ues of qq. By letting Rn be the outer radius of the first shell and plotting the 
ratio Ru/Rq in Figure[9l we see that in the limit 2/3(0) 0, Ru/Rq — > 1, 
i.e., we find numerical support that the Einstein- Vlasov-Maxwell system 
admits arbitrarily thin shell solutions. For values of go ^ 1) we see that 
F (2/3(0)) is monotonically decreasing, and as 2/3(0) approaches unity, F ap- 
proaches zero. This, however, is not the case for qo > 1 where F(2/3(0)) will at 
some point start increasing. For values of qo slightly larger than 1, F will not 
increase rapidly enough for (j5.5|) to once again be saturated. For larger val- 
ues of qo, ()5.5p will however be saturated a second time, cf. Figure fTOl where 
the graph for the case with qo = 1.25 is depicted using a different scale on 
the F— axis. From Figure [TT] it is clear that this occurs when Q = M = R. 
Hence, in the charged case the class of saturating solutions is not unique. 
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Figure 7: F as a function of 2/3(0) for a shell configuration 



k=0, 1=1, L|j=0 
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Figure 8: F as a function of 2/3(0) for a ball configuration 
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q„=1.25, k=1, 1=1, L^=10 




Figure 9: Ru/Rq as a function of ^3(0) 



Figure [12] displays the graph of p for a solution which nearly saturates the 
inequality and such that M, Q and R are almost equal, and we see that it 
is indeed not a thin shell solution. In equation (j5.2|) the quantities nii and 
ruq were defined, and roughly they represent the parts of the gravitational 
mass induced by p and q respectively. However, it should be noted that 
the nonlinearity of the equations make it impossible to completely separate 
the influences of these quantities. In figure [T2] these quantities are plotted 
in three different cases. Note in particular that in the third case, which 
corresponds to the saturating solution for which M = Q = R, rrii and mq 
are almost equal at the outer boundary of the solution. 

Figure [8] shows r(y3(0)) for the family of ball configurations with {k = 
0, / = 1, Lo = 0), for qa = {0, 0.25, 0.50, 0.75, 0.999, 1.001}. 

Although no longer monotonically decreasing, r(y3(0)) approaches zero 
as 2/3(0) approaches unity for go < 1) as in the case for the above shell 
configurations. 

6 Spirals in the mass radius diagram 

In this section we study the behavior of the total gravitational mass M and 
outer radius of the support R for one parameter families of steady states. 
The parameter is 2/3(0) while k,l,LQ and go are kept constant. In [5] it has 
been shown that in the isotropic (i.e., / = Lq = 0) and chargeless case, 
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Figure 11: Q/M and M/R as functions of y^{<d) 
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Figure 12: The profile of a saturating solution which is not a thin shell 
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Figure 13: The relation between nii and niq 
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Figure 14: {R, M)-spiral, qo = 

(R, M) forms a spiral. Furthermore, numerical evidence is given that this 
is also the case for non-isotropic solutions. Charged matter, as studied in 
this paper, displays the same behavior and changing qq merely deforms the 
spirals. This can be seen in Figures [T^ and [TSl showing the {R, M)-spirals for 
qo = and qo = 0.5 with (k = 1,1 = 5, Lo = 2). We see that increasing qo 
from to 0.5 does not change the shape significantly. These characteristics 
are displayed for all combinations of k,l and Lq- The sharp corner in these 
mass-radius spirals is a genuine feature, since for Lq > the radius of the 
support as 1/3(0) varies can change discontinuosly due to new shells which 
appear, cf. [5i p. 1829]. 

So far in this paper we have only used the ansatz (j3.14p for the distri- 
bution function. Only the analogous ansatz is used in ^5], so it might be 
interesting to see if the {R, M)-spirals are also a feature of other ansatz func- 
tions of the form f{r,w,L) = (p{E/Eo,L). For instance, with (p{E/ Eo, L) = 
sin^ {k{l - E/Eo)+) sin^ {1{L - Lo)+), and (A: = 2, Z = 3, Lq = 1, qo = 0.5), 
we get the spiral in Figure [1^1 For all cases tested, we do in fact get 

M)-spirals. We conclude that (i?, M)-spirals are a general feature of the 
Einstein- Vlasov-Maxwell system with the ansatz f{r,w,L) = (j){E/Eo,L). 
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